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Abstract 

The  operating  conditions  of  a  direct  methanol  fuel  cell  (DMFC)  are  important  with  regard  to  the  power  output  and  the  efficiency  of  the 
cell.  Especially  the  methanol  concentration  influences  the  overpotentials  of  both  the  anode  and  the  cathode.  Furthermore,  the  efficiency 
depends  on  the  amount  of  methanol  losses  due  to  methanol  permeation.  The  goal  of  this  paper  is  to  investigate  the  anode,  the  cathode  and 
the  methanol  permeation  of  a  DMFC  experimentally  to  obtain  a  set  of  parameters  for  a  DMFC  model.  The  developed  model  predicts  the 
U/I  characteristics  as  well  as  the  methanol  permeation  for  different  operating  conditions. 
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1.  Introduction 

Direct  methanol  fuel  cells  (DMFC)  are  expected  to  ap¬ 
pear  in  the  market  within  the  next  few  years  especially  for 
portable  electronic  devices  and  other  applications  in  the  low 
power  range.  There,  they  offer  the  advantage  of  easy  storage 
and  easy  refilling  of  the  liquid  methanol.  Recently,  DMFC 
stacks  and  system  in  the  power  range  from  several  watts  to 
2kW  class  have  been  reported  [1-3]. 

As  compared  to  hydrogen  driven  systems,  the  system  de¬ 
sign  and  the  operation  of  DMFC  is  much  simpler.  Although 
H2-PEFC  and  DMFC  have  nearly  the  same  theoretical  open 
circuit  voltage  (1 .23  V  versus  1 .21  V)  the  cell  voltage  and  the 
power  density  in  DMFCs  are  lower  due  to  electrochemical 
losses.  At  present,  catalyst  loadings  of  the  DMFC  electrodes 
are  substantially  higher  than  those  in  PEFC.  Therefore,  a 
present  focus  in  research  is  the  design  of  DMFC  stacks  and 
membrane  electrode  assemblies  (MEAs)  with  reduced  cat¬ 
alyst  loadings.  In  addition,  the  knowledge  about  the  pro¬ 
cesses  in  the  DMFC  as  anode  and  cathode  overpotentials  or 
methanol  permeation  is  important,  as  the  different  processes 
interact  and  may  affect  the  power  and  the  efficiency  of  the 
stacks  [4,5]. 

In  the  low  temperature  range,  DMFCs  are  operated  with 
a  liquid  water/methanol  mixture.  The  cathode  is  fed  with 
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oxygen  or  air.  In  a  direct  methanol  fuel  cell,  the  following 
catalytically  activated  reactions  take  place: 

anode  :  CH3OH  +  H20  =>■  C02  +  6H+  +  66- 

cathode  :  6H+  +  6e_  +  1.502  =>  3H20 

overall  reaction  :  CH3OH  +  1.502  =+  CO2  +  2H2O 

The  main  problem  of  the  DMFC  anode  is  the  multistage 
nature  of  the  reaction.  In  DMFCs,  CO  or  COH  occurs  as 
stable  and  adsorbed  intermediates  of  methanol  oxidation. 
These  adsorbants  lead  to  considerable  high  anodic  over¬ 
potentials  lowering  the  cell  voltages  significantly  below 
the  theoretically  expected  values.  At  anodic  potentials  less 
than  +500  mV  versus  the  reversible  hydrogen  potential, 
the  oxidation  of  methanol  is  initiated  by  a  dehydrogena¬ 
tion  reaction  of  the  methanol  [6-14].  At  higher  potentials, 
methanol  oxidation  takes  place  through  a  reaction  with  ad¬ 
sorbed  oxygen  or  OH  [6-14].  In  liquid  operated  DMFC, 
the  methanol  concentration  has  specific  importance  with 
regard  to  the  diffusion  overpotentials.  Too  low  a  methanol 
concentration  leads  to  increased  diffusion  overpotentials. 
Generally,  diffusion  overpotentials  can  be  prevented  by  the 
use  of  sufficiently  high  methanol  concentrations.  On  the 
other  hand,  this  may  cause  high  methanol  permeation  rates 
and,  therefore,  methanol  losses.  If  the  cathode  catalyst  also 
promotes  the  methanol  oxidation,  the  formation  of  a  mixed 
methanol/oxygen  potential  lowers  the  cathode  potential. 
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Nomenclature 

a 

constant 

c 

concentration 

£‘ref 

reference  concentration 

d 

membrane  thickness 

dvef 

reference  membrane  thickness 

DMFC 

direct  methanol  fuel  cell 

E 

voltage 

£o 

open  circuit  voltage 

F 

Faraday’s  constant 

fl-4 

fitting  parameters 

i 

current  density 

exchange  current  density 

fion 

ionic  current  density 

/ 

current 

K 

constant 

4nem 

membrane  thickness 

lea 

thickness  of  catalyst  layer 

m,2 

fitting  parameters 

M 

molar  concentration 

MEA 

membrane  electrode  assembly 

^clrag 

electro-osmotic  drag  factor 

PEM 

polymer  electrolyte  membrane 

R 

molar  gas  constant 

REC 

Reference  Electrode  (cathode) 

T 

temperature 

a 

transfer  coefficient 

y 

reaction  order 

7?  an 

anode  overvoltage 

0c 

cathode  overvoltage 

7?  ohm 

ohmic  losses 

®  mem 

membrane  conductivity 

The  goal  for  the  model  is  to  study  optimized  operating 
conditions  using  reliable  information  about  the  processes  in 
DMFC  as  overpotentials  or  methanol  permeation  obtained 
in  measurements.  In  addition,  the  model  can  be  used  to  ex¬ 
amine  possible  aspects  like  membrane  materials  with  differ¬ 
ent  proton  conducting  behaviour  or  different  methanol  per¬ 
meation  properties.  In  the  literature,  many  models  describ¬ 
ing  the  DMFC  are  present.  In  [15,16],  the  main  effects  in  a 
DMFC  are  described  by  mathematical  models  for  the  mass 
transport  in  the  porous  electrode  structures  and  the  poten¬ 
tial  and  concentration  distributions  in  the  electrode  regions. 
There,  the  voltage  loss  due  to  methanol  permeation  is  as¬ 
sumed  to  be  proportional  to  the  methanol  permeation.  In 
another  mathematical  approach  [17],  special  focus  is  taken 
on  the  mass  transport  of  methanol  in  the  gas  diffusion  layer 
and  the  ionic  phase  to  predict  the  methanol  permeation.  Fur¬ 
thermore,  the  mixed  potential  at  the  cathode  is  calculated  as 
a  superposition  of  the  methanol  oxidation  and  the  oxygen 
reduction.  A  macrohomogeneous  porous  electrode  describ¬ 
ing  the  reactions  and  the  transport  within  the  catalyst  zone 


of  the  methanol  electrode  is  shown  in  [18].  There,  effects  of 
catalyst  zone  thickness,  catalyst  loading  and  polymer  elec¬ 
trolyte  loading  are  discussed  by  the  means  of  a  mathematical 
approach.  In  [19],  an  equation  is  given  which  is  based  on 
a  semi-empirical  approach.  Methanol  oxidation  and  oxygen 
reduction  kinetics  are  combined  with  effective  mass  trans¬ 
port  coefficients  for  the  fuel  cell  electrodes.  Reference  [20] 
describes  a  model  particularly  for  the  DMFC  anode  tak¬ 
ing  into  account  non-Tafel  kinetics  of  electrochemical  reac¬ 
tion  of  methanol  oxidation,  diffusion  transport  of  methanol 
through  the  backing  layer  and  methanol  permeation  where 
the  total  limiting  current  density  appears  to  be  a  combina¬ 
tion  of  reaction-  and  diffusion-limiting  current  densities.  An¬ 
other  mathematical  model  for  the  anode  of  a  direct  methanol 
fuel  cell  is  presented  in  [21].  This  model  considers  the  mass 
transport  in  the  whole  anode  compartment  and  the  proton 
exchange  membrane  (PEM),  together  with  the  kinetic  and 
ohmic  resistance  effects  through  the  catalyst  layer.  Further¬ 
more,  the  influence  of,  e.g.  catalyst  layer  parameters  on  the 
methanol  permeation  and  anode  performance  is  investigated. 

2.  Experimental 

The  following  characteristics  of  DMFC  fuel  cell  have  been 
experimentally  determined: 

•  current  versus  anode  potential, 

•  cell  potential  versus  current, 

•  methanol  permeation. 

For  the  measurement  of  the  different  characteristics,  dif¬ 
ferent  test  rigs  were  used. 

2.7.  Test  cell  for  DMFC  anode  characterization 

For  the  evaluation  of  DMFC  anode,  a  separate  appara¬ 
tus  with  a  three  electrode  setup  was  used  [22].  Schematic 
drawings  of  the  measuring  apparatus  and  the  MEA  are 
shown  in  Fig.  1.  All  the  electrochemical  measurements, 
i.e.  quasi-stationary  current/potential  curves  (scan  rate: 
0.5  mV/s),  cyclic  voltammograms  (scan  rate:  20mV/s) 
and  impedance  spectra  were  performed  with  a  three  elec¬ 
trode  arrangement.  The  impedance  spectra  were  used  for 
IR- correction  and  will  not  be  discussed  in  detail.  The  elec¬ 
trode  potential,  U  (REc),  is  measured  versus  the  reference 
potential  of  the  hydrogen  reaction  established  at  the  ring 
(and  the  counter  electrode). The  anode  chamber  was  always 
supplied  with  an  aqueous  solution  of  methanol  and  the 
cathode  chamber  was  filled  with  water,  which  was  saturated 
with  hydrogen.  All  measurements  were  carried  out  under 
ambient  pressure.  Both  liquids  were  pumped  from  two  ther- 
mostated  vessels  to  the  electrode  chambers  and  backwards 
with  a  rate  of  200  ml/min.  The  temperature  was  varied  from 
30  to  70  °C.  All  the  electrodes,  which  consisted  of  both  a 
catalyst  and  a  backing  layer,  were  electrically  contacted  by 
platinum  grids  and  wires.  The  Pt  grids  were  pressed  onto 
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Fig.  1.  Three  electrode  measuring  setup  for  half  cell  measurements.  Top: 
test  equipment.  Bottom:  three  electrode  membrane  electrode  assembly  for 
half  cell  measurements. 


the  electrodes  by  using  highly  porous  filter  plates  made  of 
glass.  The  Nafion  117  membrane  had  a  diameter  of  80  mm 
and  the  electrodes  had  a  diameter  of  15  mm  corresponding 
to  an  electrode  area  of  about  1.8  cm2.  The  ring-shaped  ref¬ 
erence  electrode  placed  on  the  cathode  side  of  the  MEA  had 
an  inner  diameter  of  20  mm  and  an  outer  diameter  of  26  mm. 

2.2.  Test  cell  and  test  rig  for  cell  characterization 

The  cells  and  membrane  electrode  assemblies  were  char¬ 
acterized  by  different  test  rigs.  For  the  evaluation  of  cur¬ 
rent/density  characteristics  as  well  as  the  measurement  of 
the  methanol  permeation  the  test  cell  and  the  test  rig,  shown 
in  Figs.  2  and  3,  were  used  [19].  The  active  area  of  the  MEA 
is  20  cm2 .  Electronic  contacting  of  the  MEA  as  well  as  the 
provision  of  the  reaction  media  is  always  effected  with  the 
aid  of  a  grid  structure.  The  cell  consists  of  titanium.  The 
non-conducting  intermediate  frame  acts  on  the  one  hand  as 
an  additional  electric  insulation  of  anode  and  cathode  and, 
on  the  other  hand,  as  a  centring  aid  during  installation  of 
the  membrane  electrode  assembly.  The  inactive  rim  of  the 
MEA  along  with  two  PTFE  sealings  serves  for  the  gas-tight 
separation  of  the  anode  and  cathode  compartments. 

The  test  rig  consists  of  an  anode  loop  for  the  methanol/ 
water  mixture  with  separation  and  removal  of  the  CO2  gas 
bubbles.  Any  methanol  reaching  the  cathode  is  converted 
into  CO2  at  least  in  the  catalytic  burner  placed  behind  the 
cathode  cell  outlet.  Then,  the  amount  of  CO2  resp.  the 
methanol  permeation  is  measured  by  an  infrared  CO2  de¬ 
tector  (Vaisala  GMM  12  A)  which  is  placed  in  the  cathode 
exhaust  gas  stream.  Depending  on  the  membrane,  the  cata¬ 
lyst  layer  structures  and  the  operating  conditions  the  amount 
of  CO2  which  passes  from  the  anode  to  the  cathode  must 
be  measured  separately.  Details  of  this  method  are  given  in 
[23]. 

The  test  rig  shown  above  can  be  used  further  for  the  in¬ 
vestigation  of  the  anode  potential  at  elevated  pressures  and 
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Fig.  2.  Test  cell  for  U/I  characteristics. 
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Fig.  3.  Schematic  of  the  test  rig  used  for  the  measurement  of  methanol  permeation  and  U/l  characteristics. 


temperatures  which  are  not  suitable  for  the  test  rig  specifi¬ 
cally  designed  for  anode  potential  measurements.  For  these 
measurements,  the  cathode  chamber  is  filled  with  hydrogen. 
Then,  additional  hydrogen  is  produced  at  the  cathode  by 
applying  a  current  through  the  cell  by  means  of  an  exter¬ 
nal  current  source.  The  method  used  is  based  on  the  hydro¬ 
gen  production  in  the  cathode  chamber  of  the  test  cell.  This 
changes  the  cathode  of  the  fuel  cell  into  a  dynamic  hydrogen 
electrode  with  an  overpotential  negligible  compared  to  the 
methanol  electrode.  As  a  result,  in  a  first  approximation  the 
resulting  cell  voltage  is  the  sum  of  ohmic  losses  in  the  mem¬ 
brane  and  the  overpotentials  in  the  methanol  oxidation  elec¬ 
trode.  The  corresponding  ohmic  losses  were  determined  by 
the  current  interruption  technique  [24] .  At  elevated  temper¬ 
atures,  the  test  rig  is  pressurized  to  keep  the  water/methanol 
mixture  liquid. 

During  the  experiments,  different  types  of  membrane  elec¬ 
trode  assemblies  have  been  used  to  investigate  if  the  model 
to  be  developed  is  suitable  for  a  wide  range  of  operating 
conditions  independent  on  the  different  MEA  preparation 
methods.  In  addition,  data  taken  from  literature  are  used  for 
model  evaluation  and  verification. 

2.3.  Membrane  electrode  assembly  type  A 

These  MEA  was  prepared  in-house  by  a  spraying  method. 
The  anode  backing  layer  consisted  of  a  microlayer  com¬ 
posed  of  uncatalysed  XC-72  active  carbon  (6  mg/cm2)  and 
15  wt.%  PTFE  (Dyneon).  The  microlayer  was  supported  on 
a  carbon  cloth  (E-TEK).  The  anode  catalyst  layer  contained 
4  mg/cm2  of  Pt/Ru  catalyst  (50:50)  on  carbon  XC-72  and 
20  wt.%  Nation  from  a  5%  Nation  solution.  The  cathode  was 
similarly  constructed  with  4  mg/cm2  of  Pt  black  with  10% 
Nation.  The  backing  layer  was  composed  of  a  hydropho- 
bized  carbon  cloth  (E-TEK,  8  mg/cm2)  and  a  8  mg/cm2  mi¬ 
crolayer  consisting  of  XC-72  and  40  wt.%  PTFE.  The  MEA 


was  fabricated  by  hot  pressing  both  electrodes  with  a  Nation 
117  membrane  at  130  °C. 

2.4.  Membrane  electrode  assembly  type  B 

The  membrane  electrode  assemblies  type  B  are  used  for 
half  cell  measurements.  The  anode  catalyst  layers  were 
prepared  by  spraying  a  mixture  of  unsupported  Pt/Ru  (1:1) 
black  (Johnson  Matthey),  an  appropriate  amount  of  5% 
Nation  solution  (Aldrich)  and  isopropanol  onto  a  backing 
layer  made  of  87  wt.%  carbon  black  (XC-72)  and  13  wt.% 
PTFE  supported  by  carbon  cloth  (E-Tek,  Inc.,  ‘A’  cloth). 
The  catalyst  loading  was  1.9 mg/cm2.  The  anodes  were  hot 
pressed  (T  =  130  °C)  onto  Nation  1 17  membranes  and  acted 
as  working  electrodes.  Pt  black  catalyst  layers  with  9  wt.% 
Nation  were  used  as  counter  and  reference  electrodes. 

2.5.  Cell  voltage 

An  example  of  U/I  characteristics  of  a  DMFC  is  shown 
in  Fig.  4.  High  methanol  concentrations  prevent  the  anode 
from  the  formation  of  significant  overpotentials.  The  max¬ 
imum  current  density  achievable  with  a  0.25  M  concentra¬ 
tion  is  200mA/cm2.  Doubling  the  concentration  to  0.5  M 
doubles  the  limiting  current.  On  the  other  hand,  attention 
must  be  paid  on  the  methanol  permeation  which — for  a 
given  temperature — depends  on  the  current  density  as  well 
as  on  the  concentration.  Using  a  2  M  methanol  concentration 
lowers  the  cell  voltages  compared  to  the  1  M  concentration 
which  is  caused  by  the  impact  of  the  methanol  permeation 
on  the  cathode  potential.  Therefore,  the  optimum  methanol 
concentration  depends  on  the  current  density.  In  the  low  cur¬ 
rent  regime,  (/  <  lOOmA/cm2)  the  0.25  M  concentration 
is  best,  whereas  at  higher  currents,  (100mA/cm2  <  i  < 
250mA/cm2),  use  of  a  0.5  M  methanol  solution  is  accom¬ 
panied  with  the  highest  power  output.If  the  current  density 
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Fig.  4.  U/l  characteristics  of  a  DMFC  at  different  methanol  concentrations; 
operating  conditions:  110°C,  cathode  and  anode  pressure  3  bars,  oxidant 

o2. 

is  further  increased,  the  optimum  methanol  concentration  is 
1M. 

In  order  to  distinguish  the  different  influences  of  the 
methanol  concentration  on  the  anode  and  the  cathode,  the 
anode  overpotential  and  the  methanol  permeation  were 
determined  separately. 

The  anode  overpotential/current  characteristics  at  110°C 
and  different  methanol  concentrations  are  given  in  Fig.  5. 
With  a  methanol  concentration  of  0.5  M,  a  limiting  current 
of  about  400  mA/cm2  can  be  observed  leading  to  a  signif¬ 
icantly  increase  in  the  anodic  overpotential  600  mV  which 
is  approximately  270  mV  higher  than  the  corresponding  an¬ 
odic  overpotential  for  a  methanol  concentration  of  1  M.  This 
is  in  agreement  with  the  behaviour  of  the  whole  cell  (see 
Fig.  4).  According  to  Kauranen  et  al.  [14],  the  limiting  cur¬ 


current  density  [  A/cm2] 


Fig.  5.  Anode  potential  as  a  function  of  concentration  at  110°C;  cathode 
and  anode  pressure  3  bars;  MEA  type  A. 


rent  of  methanol  oxidation  on  platinum  ruthenium  can  pre¬ 
dominantly  be  attributed  to  a  saturation  of  the  coverage  by 
hydroxyl  groups  on  the  catalyst  surface.  However,  perme¬ 
ation  experiments  prove,  that  also  mass  transport  limitations, 
i.e.  the  diffusion  of  methanol  plays  an  important  role  (see 
below). 

In  1  and  1.5  M  methanol  solution,  the  limiting  current  is 
much  higher  than  the  maximum  current  density  investigated. 
Furthermore,  in  case  of  the  1.5  M  methanol  solution,  no 
reaction  or  diffusion  overpotential  is  observed. 

Finally,  the  methanol  permeation  for  the  given  operating 
conditions  is  shown  in  Fig.  6.  Three  different  methanol 
concentrations  (0.5,  1  and  1.5  M)  have  been  tested.  The 
measurements  include  corrections  for  the  CO2  passing  from 
the  anode  to  the  cathode  and  are  detected  additionally  in  the 
previously  described  CO2 -sensor.  Without  CO2 -correction, 


Fig.  6.  Methanol  permeation  as  a  function  of  methanol  concentration  and  current  density. 
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the  measured  methanol  permeation  is  too  high  especially 
at  high  current  densities.  Details  of  the  method  are  given 
in  [23].  Using  the  lowest  methanol  concentration  of  0.25  M 
the  methanol  permeation  decreases  to  approximately 
20mA/cm2  at  a  current  density  of  400mA/cm2  indicating 
a  depletion  of  methanol  in  the  interface  between  anode  cat¬ 
alyst  layer  and  membrane.  This  is  in  good  agreement  with 
the  observed  limiting  current  behaviour  at  the  same  current 
density. 

Furthermore,  the  methanol  permeation  as  a  function  of 
temperature  is  given  in  Fig.  7  for  the  open  circuit  condition. 
Increasing  the  temperature,  e.g.  from  60  to  110°C  doubles 
the  methanol  permeation.  This  is  mainly  due  to  increased 
diffusivity  of  methanol  in  the  water  and  the  membrane.  With 
rising  current  densities,  the  electro-osmotic  drag  becomes 


predominant  compared  to  the  diffusion.  The  drag  coeffi¬ 
cient  ftdrag,  describing  the  ratio  of  transported  molecules  to 
transporting  protons,  is  temperature  dependent  as  well  e.g. 
/i drag  =  2  at  15  °C  reaching  5.1  at  130  °C  [25]. 

To  investigate  the  influence  of  the  temperature  on  the 
anode  overpotential  characteristics,  measurements  in  the  test 
cell  with  three-electrode  ME  As  were  performed.  The  results 
of  these  measurements  are  shown  in  Fig.  8. 


3.  Model  development 

The  new  model  is  based  on  the  experiments  presented 
above  with  the  processes  schematically  shown  in  Fig.  9. 
The  goal  of  this  model  is  to  predict  the  main  processes  as 


Fig.  8.  IR- corrected  anode  polarization  curves  as  a  function  of  temperature  at  1  M  methanol  concentration.  MEA  type  B. 
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Fig.  9.  Schematic  of  a  direct  methanol  fuel  cell. 


methanol  permeation  and  overpotentials  resulting  in  specific 
current/voltage  characteristics  in  a  wide  range  of  operating 
conditions.  In  a  future  step,  the  model  will  be  used  to  predict 
the  power  and  efficiency  of  a  DMFC  system  as  a  function  of 
its  operation  conditions.  For  reasons  of  clarity  and  usability 
the  model  should  be  simplified  to  contain  as  few  parameters 
as  possible. 

The  cell  voltage  of  a  DMFC  is  calculated  by  the  theoretical 
open  voltage  less  the  overpotentials  q an,  ijc  and  the  ohmic 
losses  ^ohm- 

E  —  E0  ?7an  Oc  7? ohm  (1) 


species  followed  by  a  second  step  containing  the  associative 
ionization  of  adsorbed  CO  with  water  molecules. 
i*cy[a  exp  (2. 303 77 /Z?) 

i  =  -  (2) 

cm  +  Azz  exp  (2. 3 03  27//?) 

with  4  as  the  exchange  current  density  of  the  anode  [A/cm2], 
cm  the  methanol  concentration  in  the  catalyst  layer,  b  rep¬ 
resents  the  Tafel  slope,  K  and  a  are  constants. 

For  large  methanol  concentrations  cm  and  small  overpo¬ 
tentials  T)  Eq.  (2)  describes  zero-order  kinetics.  In  this  case, 
the  exponent  in  the  denominator  can  be  neglected  with  the 
result  that  the  current  i  is  independent  of  the  methanol  con¬ 
centration  and  is  a  function  only  of  the  overpotential. 

/  2.303q 
i  =  i*a  exp  I  - 

V  b 

for  large  methanol  concentrations  and  small  q  (3) 


In  contrast,  for  high  overpotentials,  q ,  the  concentration 
term  in  the  denominator  can  be  neglected  indicating  a  tran¬ 
sition  to  first-order  kinetics.  Then,  the  current  is  limited  by 
the  finite  rate  of  methanol  adsorption  on  the  catalyst  layer 
which  is  proportional  to  the  methanol  concentration. 


i*cm 

~~K~ 


limiting  current  for  high  q 


In  the  following,  both  the  physical  model  and  the 
semi-empirical  model  are  fitted  to  the  experimentally  ob¬ 
served  data.  The  fitting  parameters  are  summarized  in 
Table  1. 


The  theoretical  open  circuit  voltage  Eq  can  be  calculated 
using  Nernst’s  equation.  Its  value  is  Eq  =  1.21  V  for  stan¬ 
dard  conditions. 

3.1.  Physical  model  of  the  anode  overpotential 

The  anodic  overpotential  is  mainly  a  function  of  the  re¬ 
action  rate,  the  methanol  concentration  and  the  temperature. 
Eq.  (2)  given  by  Meyers  and  Newman  [26]  is  a  simplification 
of  the  four  step  mechanism  proposed  by  Gasteiger  et.  al  in 
[11]  resulting  in  a  simplified  two  step  mechanism.  The  first 
step  is  adsorption  of  methanol  and  ionization  of  adsorbed 


3.2.  Semi-empirical  equation 

The  anode  overpotential  is  fitted  by  the  following 
semi-empirical  equation: 

??an  =  f\+h  In (i)  +  h  exp  ( — )  [V]  (5) 

\CM/ 

This  equation  delivers  the  overvoltage  directly  as  a  func¬ 
tion  of  the  current  density  which  is  beneficial  for  numerical 
implimentation  to  calculate  Eq.  (1),  whereas  the  equation  of 
Meyers  and  Newman  cannot  be  transformed  to  give  f]an  as 
a  function  of  the  current  density  i. 


Table  1 

Fitting  parameters  for  the  polarization  curves  of  the  MEA  anodes  shown  in  Fig.  10 


Anode  fitting  parameters  for 


Concentration  (M)  The  semi-empirical  equationa  Meyers  and  Newmanb 


/l  (V) 

fl  (V) 

/3  (V) 

f\  (mol/1)/ (A/cm2) 

4  (A/cm2) 

a  [-] 

b  (V) 

K  (M) 

MEA  type  A  110°C,  3  bars 

0.5-1. 5  0.3949 

0.0603 

4.3565E-5 

10.81 

1 .22E-4 

0.2 

0.074 

1.6E-4 

MEA  described  in  literature  [14] 
0.02-0.5  0.594 

0.0369 

4.3565E-5 

8.4 

5.11E-7 

0.2 

0.085 

4.955E-7 

a  11a  =  fi  +  h  ln(0+/3  exp  El  . 

b  •  _  i*  CMflexp(2.303??/fr) 
CM+Kaexp(2.303ri/b) ' 
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In  general,  this  equation  is  a  combination  of  Tafel-kinetics 
with  an  additional  exponential  term  describing  the  formation 
of  limiting  current  behaviour. 

The  constant /i  and  the  logarithmic  expression  of  Eq.  (5) 
can  be  derived  partially  from  Eq.  (2),  respectively  Eq.  (3) 
for  the  case  of  small  overpotentials  and  large  methanol  con¬ 
centrations.  Then,  they  contain  the  tafel  slope  b  as  well  as 
the  exchange  current  density  /* .  The  logarithmic  term  dom¬ 
inates  in  the  regime  of  low  current  densities. 

Of  course,  alternatively  these  parameters  can  be  obtained 
by  simple  fit  procedures  on  the  basis  of  the  experimental 
data.  This  alternative  approach  is  necessary  if  the  first  ap¬ 
proach  by  Meyers  and  Newman  does  not  give  sets  of  param¬ 
eters  which  fit  the  experimental  data.  Comparing  the  coeffi¬ 
cients  in  Eq.  (3)  and  Eq.  (5)  leads  to  the  following  analogy: 


[V]  /i 


b 

2.303’ 


[V]  h 


b 

2.303 


In 


A/cm2 

a/* 


The  last  term  of  the  sum  in  Eq.  (5)  finally  increases  the 
overpotential  rjan  with  increasing  current  density.  It  can  be 
easily  seen  that  higher  methanol  concentrations  shift  the 
begin  of  limiting  current  behaviour  towards  higher  current 
densities. 

A  similar  expression  is  described  in  [27]  for  the  cathodic 
overpotential  of  a  PEM  with  the  exponential-term  for  fit¬ 
ting  the  limiting  current  behaviour  due  to  oxygen  depletion. 
Analogously,  the  exponential-term  is  used  in  our  DMFC 
modeling  to  describe  the  increase  of  the  anode  overpotential 
at  too  low  methanol  concentrations. 

As  experimentally  observed,  the  anode  overpotential  of  a 
DMFC  follows  typical  curves  with  characteristics  as  asymp¬ 
totes,  ranges  with  semi-linear  behaviour,  inflection  points 
and  offset  to  the  hydrogen  reversible  potential.  For  a  given 
DMFC  anode,  these  characteristic  points  mainly  depend  on 
the  temperature  and  the  methanol  concentration. 


Fig.  10  shows  experimental  data  and  fitting  curves  for  an 
in-house  prepared  MEA  (left)  and  for  data  taken  from  the 
literature  (right,  [14]).  The  fitting  parameters  are  summa¬ 
rized  in  Table  1 .  For  each  MEA  investigated,  the  parameter 
f\  is  reciprocal  to  the  concentration,  i.e.  the  impact  of  dif¬ 
fusion  and  reaction  overpotentials  is  linearly  coupled  to  the 
methanol  concentration  if  the  temperature  is  maintained.  Of 
course,  changing  the  temperature  influences  the  beginning 
of  limiting  current  behaviour  as  increasing  temperature  in¬ 
creases  the  diffusion  coefficient. 

To  investigate  the  influence  of  the  temperature  on  the 
current/density  characteristics,  the  measurements  described 
above  with  constant  methanol  concentration  and  tempera¬ 
tures  varying  in  the  range  30-70  °C  were  analysed  and  fit¬ 
ted.  As  shown  in  Fig.  11,  the  potential/current  curve  shifts 
nearly  parallel  towards  lower  overpotentials  if  the  temper¬ 
ature  increases.  This  is  due  to  increased  kinetics  and  de¬ 
creased  activation  overpotential  of  the  methanol  oxidation. 
Furthermore,  the  beginning  of  the  limiting  current  behaviour 
is  shifted  to  higher  current  densities.  The  fitting  parame¬ 
ters  regarding  the  temperature  variation  are  summarized  in 
Table  2.  Although  different,  but  similar,  ME  As  were  inves¬ 
tigated  the  fitting  parameters  f\  and  f\  show  a  clear  trend 
concerning  the  temperature  which  can  be  used  in  the  sim¬ 
ulation  studies.  In  the  temperature  between  3  and  1 10  °C, 
the  parameter /i  describing  the  offset  depends  linearly  on 
the  temperature  whereas  the  parameter  fa  shows  Arrhenius 
behaviour  (Fig.  12). 

The  corresponding  fitting  equations  are: 

(T  —  273.15AA 

fi  =  -0.0021  f - — - J  +  0.6132  (7) 


f\  =  0.0297  exp 


2261.7 


T 


Fig.  10.  /^-corrected  polarization  curves  of  different  DMFC  anodes  as  a  function  of  methanol  concentration  (dots:  experimental  data,  solid  lines:  fit  by 
the  semi-empirical  equation,  dashed  lines:  fit  by  the  equation  of  Meyers  and  Newman).  Left:  MEA  type  A  at  higher  methanol  concentrations.  Operating 
conditions:  temperature  110°C,  pressure  3  bars,  high  overstoichiometric  flow  of  methanol.  Right:  experimental  data  (dots)  taken  from  Kauranen  et  al. 
[14]  at  60  °C  in  1M  H2SO4,  catalyst  loading  7.9  mg/cm2  Pt/C,  ambient  pressure. 
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current  density  i  [A/cm2] 


Fig.  11.  /^-corrected  polarization  curves  as  a  function  of  temperature,  DMFC  MEA  type  B;  dots:  experimental  data,  solid  lines:  fit  by  the  semi-empirical 
equation,  dashed  lines:  fit  by  the  equation  of  Meyers  and  Newman. 


Table  2 


Fitting  parameters  for  different  temperatures  at  1  M  methanol  solution 


Temperature  (°C) 

Anode  fitting  parameters  for  the 

semi-empirical  equation3 

Anode  fitting  parameters  for  Meyer s-Newmanb 

fl  (V) 

fa  (V) 

fa  (V) 

h  (mol/l)/(A/cm2) 

i*  (A/cm2) 

a  [-] 

b  (V) 

K  (M) 

MEA  type  B  1  M, 

ambient  pressure 

30 

0.5679 

0.05525 

4.3565E-5 

53.69 

3.97E-6 

0.2 

0.0850 

2.15E-8 

40 

0.5286 

0.05525 

4.3565E-5 

40.62 

4.95E-6 

0.2 

0.0795 

1.90E-9 

50 

0.5027 

0.05525 

4.3565E-5 

30.41 

5.96E-6 

0.2 

0.0748 

1.60E-8 

60 

0.4785 

0.05525 

4.3565E-5 

25.29 

6.96E-6 

0.2 

0.0705 

1.45E-8 

70 

0.4629 

0.05525 

4.3565E-5 

23.51 

7.0E-6 

0.2 

0.0690 

1.30E-8 

MEA  type  A  1  M,  3  bars 

110 

0.3949 

0.0603 

4.3565E-5 

10.81 

1.22E-4 

0.2 

0.074 

1.6E-4 

a  r] a  =  fi  +  h  ln(0+/3  exp^. 

b  •  i *  cy[atxp{2.303r]/b) 

1  CM+Kaexp(2.303ri/b)  ' 


1/T  [1/K] 


Fig.  12.  Influence  of  the  temperature  on  the  fit  parameters  f\  and  fa  of  the  semi-empirical  equation  for  two  different  MEAs.  MEA  type  A  (□);  MEA 
type  B  (♦). 
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The  parameters  and  do  not  depend  on  the  tempera¬ 
ture.  For  a  given  MEA,  these  values  can  be  set  constant  for 
all  operating  conditions. 

3.3.  Evaluation  of  the  two  different  anode  fitting 
approaches 

Comparing  the  two  modeling  approaches — the  semi- 
empirical  approach  and  the  equation  of  Meyers  and 
Newman — it  can  be  seen  that  both  approaches  are  in  general 
suitable  for  modeling  the  typical  characteristics  of  the  ex¬ 
perimentally  determined  anode  overpotential.  In  particular, 
both  approaches  include  the  limiting  current  behaviour  as  a 
function  of  the  methanol  concentration  and  the  Tafel-like  ki¬ 
netics  at  low  overpotentials.  With  regard  to  the  congruence 
of  experimental  data  and  fittig  curves,  the  equation  of  Mey¬ 
ers  and  Newman  fits  best  for  low  concentrations  and  low 
temperatures.  On  the  other  hand,  for  higher  temperatures 
(>50  °C)  it  predicts  too  abrupt  limiting  current  behaviour. 
As  DMFC  are  operated  at  elevated  temperatures  to  obtain 
high  power  densities,  we  have  chosen  the  semi-empirical 
approach. 

3.4.  Modeling  the  methanol  permeation 

For  a  calculation  of  the  cathode  potential,  the  methanol 
permeation  must  be  taken  into  account.  The  methanol 
permeation  is  the  sum  of  a  diffusive  and  a  convective 
(electro-osmotic)  fraction.  In  the  following,  the  absence  of 
a  pressure  gradient  between  anode  and  cathode  is  assumed. 
For  given  operating  conditions  as  temperature  and  pressure, 
the  local  methanol  permeation  is  assumed  to  be  a  linear 
function  of  the  local  current  density  with  an  offset  m\  and 
the  slope  m2. 

Jperm  —  m\  +  m2i  (9) 

The  semi-empirical  expression  Eq.  (9)  is  a  result  of  dif¬ 
ferent  influences.  At  first,  the  anode  catalyst  layer  repre¬ 
sents  a  sink  for  methanol  and  at  second,  the  methanol/water 
mixture  at  the  boundary  between  anode  catalyst  layer  and 
membrane  is  transported  not  only  by  diffusion  but  also  by 
electro-osmosis  to  the  cathode.  Therefore,  with  increasing 
current  density,  the  methanol  concentration  in  the  bound¬ 
ary  between  catalyst  layer  and  membrane  decreases  with  the 
result  that  both  the  concentration  gradients  between  anode 
and  cathode  decrease  as  well.  Consequently,  a  decrease  in 
methanol  permeation  should  be  expected  for  increasing  cur¬ 
rent  densities. 

But,  on  the  other  hand,  increasing  the  current  density  leads 
to  an  increase  of  the  methanol  transport  by  electro-osmotic 
drag.  In  certain  conditions,  the  consumption  of  methanol 
in  the  anode  catalyst  layer  can  be  overcompensated  by  the 
electro-osmotic  transport.  Experiments  show  that  at  low  con¬ 
centrations  e.g.  below  2M  at  60  °C,  the  slope  is  negative, 
i.e.  the  methanol  permeation  decreases  with  increasing  cur¬ 
rent  density  due  to  the  methanol  consumption  in  the  catalyst 


layer.  At  concentrations  above  2M,  the  methanol  perme¬ 
ation  increases  as  well  because  the  anode  catalyst  layer  has 
a  limited  capacity  towards  methanol  oxidation  which  can  be 
exceeded  with  increasing  concentration.  Consequently,  the 
slope  m2  decreases  with  increasing  drag  coefficient. 

The  offset  m\  describes  the  methanol  permeation  at  open 
circuit  conditions.  It  mainly  depends  on  the  methanol  con¬ 
centration  gradient  between  anode  and  cathode  which  is 
maintained  by  the  continuous  oxidation  of  methanol  at  the 
cathode.  For  reasons  of  simplification,  the  concentration  is 
assumed  to  be  proportional  to  the  thickness  of  the  membrane 
without  taking  concentration  gradients  e.g.  in  the  diffusion 
layers  into  account.  As  the  concentration  gradient  in  the  dif¬ 
fusion  layer  is  known  to  be  nearly  negligible  [15,16,20,28] 
the  error  made  by  this  assumption  is  reasonable  small. 

mi  =  4.65 E  -  2  —  —  exp(0.01329(T  -  273.15^))-^ 

dY&  f  1M  cm2 

(10) 


The  thickness  dYQ f  is  the  reference  thickness  of  the  Nation 
117  membrane  used  in  the  experiments  (thikness  175  jam). 

The  drag  coefficient  ^drag  is  important  for  the  slope  m2. 
It  describes  the  transport  of  molecules  by  the  movement  of 
the  ions  towards  the  cathode. 

..drag  _  bon  .... 

7VH20+Me0H  —  ^drag  ^  UU 


The  slope  m2  represents  the  influence  of  the  current  den¬ 
sity  on  the  methanol  permeation  and  summarizes  effects  by 
methanol  consumption  in  the  anode  and  methanol  transport 
by  electro-osmosis.  It  is  a  function  of  the  methanol  con¬ 
centration  and  the  drag  coefficient  ^drag-  For  m2  <  0  the 
methanol  permeation  decreases  with  increasing  current  den¬ 
sity.  Finally,  m2  =  0  indicates  a  constant  methanol  perme¬ 
ation  for  all  current  densities.  If  the  methanol  concentration 
is  substantially  high  e.g.  higher  than  2  M,  the  methanol  per¬ 
meation  increases  with  increasing  current  density.  Accord¬ 
ingly,  high  drag  coefficients  ^drag  have  the  same  impact  as 
high  methanol  concentrations,  i.e.  high  drag  coefficients  lead 
to  increasing  methanol  permeation.  For  the  most  commonly 
used  membrane  materials  like  Nafion  and  concentrations  in 
the  range  of  1  M,  the  slope  is  negative,  i.e.  increasing  current 
density  lowers  the  methanol  permeation. 

Based  on  these  facts,  the  following  semi-empirical  ex¬ 
pression  is  assumed  for  m2 : 


m2  =  —0.4  +  0.27 


^drag  C 


"drag  Cref 


(12) 


The  value  ^drag  specifies  the  number  of  entrained 
molecules  per  proton.  The  drag  coefficient  as  function 
of  temperature  was  measured  with  different  experimen¬ 
tal  setups  [25,29],  and  is  in  the  range  of  2H20/H+  at 
15  °C-5.1  H20/H+  at  130  °C  for  Nafion.  Using  the  re¬ 
sults  in  [25]  the  linearized  expression  for  the  temperature 
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Fig.  13.  Fit  function  (solid  line)  for  the  measured  methanol  permeation  (dots)  in  dependence  of  the  current  density  at  different  methanol  concentrations. 


dependency  is 

ndiag=4+-^L(T-363K)  (13) 

Fig.  13  shows  the  methanol  permeation  at  different 
methanol  concentrations.  With  regard  to  the  0.5  and  1  M 
concentrations,  the  fit  is  reliable.  For  the  highest  methanol 
concentration,  1.5  M,  the  fit  slightly  underestimates  the  in¬ 
crease  of  methanol  permeation  at  higher  current  densities. 
For  the  practical  evaluation  of  fuel  cell  operating  conditions 
this  error  is  not  relevant  as  the  concentration  1.5  M  will  not 
be  adjusted  for  practical  use  as  it  is  accompanied  by  too 
high  methanol  losses. 


3.5.  Cathode  overpotential 


The  potential  of  the  cathode  is  influenced  by  the  methanol 
permeation,  the  temperature  as  well  as  the  current  density. 
The  expression  for  the  current  density  i  for  the  oxygen  elec¬ 
trode  without  permeation  is  derived  from  Tafel  equation 


c02 

CC>2,ref 


y  (olF  \ 

“P^J 


(14) 


Table  3 

Fitting  parameters  of  Eq.  (17)  for  the  U/I  characteristics  shown  in  Fig.  14, 
simulating  a  DMFC  at  3  bars  oxygen  and  110°C 


Cathode  fitting  parameter 


Numerical  expression 


RT 

aF 

RT 


In 


A/cm 
Cat  6 


2  1 


In 


olF  |_A/crm 


c02,ref 
cC>2 


r 


0.343  V(0.380V  at  85  °C) 


0.0061n 


i 

A/cm2 


l 


0.5  V 


V 

V 


k(i) 


0.15  V  + 


A/cm2  /  |_A/cm2_ 


For  reasons  of  simplification,  the  methanol  permeation  is 
regarded  as  an  additional  impact  which  increases  the  cathode 
overpotential  resulting  in: 


RT 

ric  =  —  In 
at 


^cat  i* 


(  cP2,ref 

V  co2 


V 


+  k(i)  ^perm 


(16) 


respectively 


Yl 


+  k(i)i 


perm 


(17) 


with  /cat  being  the  catalyst  layer  thickness,  /*  the  exchange 
current  density,  co2  the  oxygen  molar  concentration,  Cq *  the 
reference  oxygen  molar  concentration  which  is  a  parame¬ 
ter  of  the  given  reaction  and  given  catalyst,  y  the  reaction 
order  (assumed  to  be  1),  xik  the  overpotential  and  a  being 
the  transfer  coefficient  [20,27].  An  important  assumption  of 
Eq.  (14)  is  that  the  catalyst  layer  is  thin  enough  to  ensure 
constant  reaction  rate,  oxygen  concentration  and  overpoten¬ 
tial.  Eq.  (14)  allows  to  express  the  overpotential  r]c  by  i: 


0c  = 


aF 


i 


f  cP2,ref 

V  co2 


Yl 


(15) 


with  k  as  an  empirical  additional  factor  which  considers  the 
impact  of  methanol  permeation  in  comparison  to  the  impact 
of  the  current  density  i  on  the  cathode  potential.  An  impor¬ 
tant  observation  in  our  modeling  studies  is  that  the  impact  of 
methanol  permeation  increases  with  increasing  current  den¬ 
sity.  As  a  result,  the  factor  k  itself  depends  on  the  current 
density.  The  factor  k  was  determined  experimentally  by  fit¬ 
ting  the  experimental  data  in  Fig.  4  resulting  in  the  empiric 
linear  expression  for  k  given  in  Table  3.  For  the  given  set 
of  U/I  characteristics  at  different  methanol  concentrations  k 
starts  at  a  relatively  low  value  of  0. 15  V  per  A/cm2  methanol 
permeation  rate.  Increasing  the  electrical  current  increases 
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k  to  e.g.  0.4  V  per  mA/cm2  methanol  permeation  rate  at  an 
electrical  current  density  of  0.5A/cm2.  For  different  tem¬ 
peratures,  the  exchange  current  density  /*  changes  as  well: 


i*(T)  =  if  exp 


(18) 


Increasing  the  temperature  from  85  to  110°C,  it  tripples 
the  exchange  current  density  which  is  a  reasonable  value  for 
this  temperature  regime.  It  corresponds  to  an  activation  en¬ 
ergy  of  50  kJ/mol.  Detailed  data  for  the  activation  energy  in 
DMFC-cathodes  are  given  in  [30].  For  reasons  of  simplifica¬ 
tion,  the  other  parameters  given  in  Eq.  (17)  are  set  constant, 
in  particular  the  temperature  influence  of  the  methanol  per¬ 
meation  on  the  cathode  overpotential  is  not  yet  considered 
due  to  the  lack  of  experimental  data. 

Finally,  ohmic  losses  are  mainly  caused  by  the  membrane 
resistance.  For  reasons  of  simplification,  other  ohmic  losses 
as  losses  in  the  gas  diffusion  layers  and  contact  resistances 
are  not  taken  into  account.  The  losses  in  the  membrane  are 
proportional  to  the  thickness  /mem- 


7?  ohm  — 


Oncm  . 
- 1 

^mem 


(19) 


The  thickness  of  the  Nation  117  membrane  is  175  pan 
with  a  bulk  conductivity  of  0.1  S/cm. 


3.6.  Model  verification  for  U/I  characteristics 


The  cathode  fitting  parameters  are  given  in  Table  3.  Us¬ 
ing  Eq.  (18)  it  can  be  seen  that  the  cathode  exchange  current 
density  at  85  °C  decreases  to  a  third  compared  to  1 10  °C.  To 
evaluate  the  model,  the  experimentally  obtained  U/I  charac¬ 
teristics  and  methanol  permeation  rates  are  compared  to  the 
simulated  data.  Fig.  14  compares  the  experimental  data  with 


Fig.  15.  U/I  characteristics  of  a  DMFC  single  cell  for  a  temperature 
variation.  Experimental  data:  dots,  fit:  lines.  MEA  type  A,  anode  and 
cathode  pressure  3  bars. 


the  simulated  data  for  a  DMFC  at  1 10  °C  and  3  bars  operat¬ 
ing  pressure  with  different  methanol  concentrations.  Above 
current  densities  of  30  mA/cm2,  the  fit  is  reliable.  The  effect 
of  methanol  permeation  is  predicted  correctly.  Increasing  the 
methanol  concentration  from  1  to  2  M  at  400  mA/cm2  low¬ 
ers  the  cell  voltage  from  400  to  300  mV.  At  50  mA/cm2,  the 
cell  voltage  decreases  only  by  50  mV  using  2M  instead  of 
1  M  methanol  solution. 

Finally,  the  influence  of  the  temperature  on  the  cell  voltage 
characteristics  is  shown  in  Fig.  15.  The  current/voltage  char¬ 
acteristics  for  85  °C  has  been  measured  at  a  methanol  con¬ 
centration  of  1 M.  The  temperature  dependent  fit  parameters, 
/i  and/4,  for  the  anode  have  been  derived  from  the  temper- 


Fig.  14.  U/I  characteristics  of  a  DMFC  single  cell.  Experimental  data:  dots  with  thin  lines,  fit:  thick  lines.  The  fit  is  reliable  for  current  densities 
>30 mA/cm2.  Fit  parameters  are:  110°C,  MEA  type  A,  anode  and  cathode  pressure  3 bars,  anode  set  of  parameters  according  to  Table  1,  cathode  fitting 
parameters  according  to  Table  3. 
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ature  dependency  given  in  Eq.  (7)  and  Eq.  (8)  Their  values 
at  85  °C  are  f\  =  0.4347  V,  /4  =  16.4  (mol/1) /(A/cm2), 
the  parameters /2  and/3  are  unchanged. 


4.  Conclusions 

A  simulation  model  has  been  developed  to  predict  the 
current/density  characteristics  of  a  DMFC  for  a  wide  range 
of  operating  conditions.  This  model  is  based  on  several 
separate  experiments  for  the  evaluation  of  the  anode  and 
cathode  overpotentials,  the  methanol  permeation  and  the 
U/I  characteristics  to  obtain  a  broad  set  of  data  describing 
the  characteristics  of  the  single  parts  of  the  membrane  elec¬ 
trode  assembly.  The  methanol  permeation  has  a  substantial 
impact  on  the  current  density  and  the  efficiency.  With  re¬ 
gard  to  the  methanol  permeation  the  model  predicts  the 
impact  of  the  mixed  potential  at  the  cathode  lowering  the 
U/I  characteristics  of  the  DMFC.  This  impact  is  taken  into 
account  by  an  empirical  factor  which  has  been  found  to 
increase  with  increasing  current  density. 

In  further  studies,  this  model  will  be  used  to  determine 
the  influence  of  design  parameters  as  membrane  thickness, 
drag  factor  or  electrical  conductivity  on  the  performance  of 
DMFC. 
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